The shape of the UV ionizing background is reconstructed 
from optically thin metal absorption-line systems identified in 
spectra of HE 2347-4342, Q 1157+3143, and HS 1700+6416 
in the redshift interval 1.8 < z < 2.9. The systems are an- 
alyzed by means of the Monte Carlo Inversion method com- 
pleted with the spectral shape recovering procedure. The UVB 
spectral shape fluctuates at 2.4 < z < 2.9 mostly due to ra- 
diative transfer processes in the clumpy IGM. At z < 1.8, 
the IGM becomes almost transparent both in the H I and He II 
Lyman continua and the variability of the spectral shape comes 
from diversity of spectral indices describing the QSO/AGN in- 
trinsic radiation. At z > 2.4, the recovered spectral shapes 
show intensity depression between 3 and 4 Ryd due to He II 

■ Lya absorption in the IGM clouds (line blanketing) and con- 
tinuous medium (true Gunn-Petersen effect). The mean He II 

< ^ • Lya opacity estimated from the depth of this depression cor- 
responds within l-2er to the values directly measured from the 

■ H i/He II Lya forest towards the quasars studied. The observed 
' scatter in rj = N(He ll)/iV(H I) and anti-correlation between 

7V(H I) and ?/ can be explained by the combined action of vari- 
^•3" able spectral softness and differences in the mean gas density 
,— ' ■ between the absorbing clouds. Neither of the recovered spectral 
shapes show features which can be attributed to the putative in- 
put of radiation from soft sources like starburst galaxies. 
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1. Introduction 

This paper continues our study of the UV ionizing background 
(UVB) of the metagalactic radiation in the range 1-10 Ryd as 
reconstructed from optically thin QSO metal absorption sys- 
tems. In the previous papers, our analysis of the absorbers at 
z w 3 (Agafonova et al., 2005, hereafter Paper I) and in the 
range 1.46 < z < 1.8 (Reimers et al., 2006, hereafter Paper II) 
has led to the following conclusions: 

1 . At z « 3, the spectral shape of the UV ionizing radiation 
shows a sharp reduction in flux in the energy range 3 Ryd 
< E < 4 Ryd which can be interpreted as manifestation 
of the He II Gunn-Petersen effect — He II Lya continuous 
absorption in the intergalactic medium. 

2. At z < 1.8, the ionizing spectra turn out to be much harder 
at E > 3 Ryd as compared to model spectra of Haardt 
& Madau (1996, hereafter HM). This indicates that inter- 
galactic clouds with jV(Hi) > 10 15 cm~ 2 responsible for 
the He II break (absorption in the He II Lyman continuum) 
become more rare at low z than was supposed by the cloud 
statistics in HM. 

3. The ionizing background both at z « 3 and at z < 1.8 is 
dominated by radiation from QSO/AGNs. The input of the 
soft radiation from galaxies/stars is negligible, limiting the 
escape fraction of the galactic UV photons to / osc < 0.05. 
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In the present paper, we reconstruct the underlying UVB using 
metal absorbers from the range 1.8 < z < 2.94 identified in 
spectra of three bright QSOs: HE 2347-4342 (z om = 2.88), Q 
1157+3143 (z cm = 2.96), and HS 1700+6416 (z cm = 2.74). In 
the considered redshift range, the shape of the UVB is closely 
related to the process of He II reionization in course of cos- 
mic time. Recent observations of the H I/He II Lya forest to- 
wards HE 2347^1342 (Shull et al. 2004; Zheng et al. 2004, 
hereafter Z04) and HS 1700+6416 (Fechner et al. 2006a, here- 
after F06) revealed (i) significant fluctuations of the ratio r\ = 
N(He n)/N(H I), (ii) anti-correlation between rj and iV(H I), 
and (Hi) decrease of median r\ towards lower z. Although a 
large part of the scatter of rj seems to be artificial and caused 
by noise in Hell data, errors in Hi and Hell column densi- 
ties, unidentified blends etc. (see F06), similarity of the re- 
sults obtained both for the HE 2347^1342 and HS 1700+6416 
sightlines supposes that spatial fluctuations of 77 do really exist. 
These fluctuations suggest a variable softness S = J912/ ' J22& 
of the metagalactic radiation field with higher r\ value corre- 
sponding to softer UVB. A diversity of the spectral indices de- 
scribing QSO continuum radiation, an input of radiation from 
soft sources of type starburst galaxies, and radiative transfer 
through clumpy media are the reasons usually considered for 
the UVB variability. Calculations of the UVB spectra account- 
ing for all these effects are still well beyond modern compu- 
tational facilities, although some attempts to explain the ob- 
served scatter of rj via the modeling of the H I/He II Lya forest 
have been made (Maselli & Ferrara 2005; Bolton et al. 2006). 
Thus, direct reconstructions of the UVB spectral shape from 
metal absorption-line systems at different redshifts and along 
different sightlines can help to clarify the character of the UVB 
fluctuations and to determine the sources of the ionizing back- 
ground. 
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The paper is organized as follows. Sect. 2 contains some 
qualitative information outlining main steps in the UVB recov- 
ering procedure. Selected absorption systems are analyzed in 
Sect. 3. The obtained results are discussed and summarized in 
Sects. 4 and 5. 

2. Recovering the shape of the underlying 
radiation field 

Computational procedure used to recover the spectral shape of 
the ionizing radiation from optically thin metal systems is de- 
scribed in Paper I. Here we briefly outline the basics needed to 
understand the results presented below in Sect. 3. 

Absorption systems are analyzed by means of the Monte 
Carlo Inversion (MCI) procedure (Levshakov, Agafonova & 
Kegel 2000, hereafter LAK; Levshakov et al. 2002, 2003a,b). 
The MCI is based on the assumption that all lines observed 
in the absorption system are formed in a cloud where the gas 
density and velocity are random functions of the coordinate x 
along the cross line. It is also assumed that the metal abundance 
throughout the cloud is constant, the gas is in thermal and ion- 
ization equilibrium and optically thin for the ionizing UV ra- 
diation. This implies that ion fractions are determined entirely 
by the gas density Ub_(x) [or, equivalently, by the ionization 
parameter U{x) oc 1/ «h] and by the spectral energy distribu- 
tion of the background ionizing radiation in the range 1-10 Ryd 
(defined by the ionization edges of the observed metal ions — 
from Si II to O VI, see Fig. 1). 

Since the gas density varies from point to point along the 
sightline, the ion fractions vary as well. Thus, for a given point 
within the line profile the observed intensity results from a mix- 
ture (superposition) of different ionization states of the gas due 
to irregular Doppler shifts of the local absorption coefficient 
(see Fig. 1 in LAK). The dependence of the ion fraction on the 
gas density (ionization curve) is different for each ion and this 
is why the line profiles differ in spite of being produced by the 
same gas cloud. 

The following physical parameters are directly estimated 
in the computational procedure: the mean ionization parameter 
/Jo, the total hydrogen column density A^h, the line-of-sight 
velocity dispersion, er v , and density dispersion, er y , of the bulk 
material [y = n-n{x) / uq\, and the chemical abundances Z a of 
all elements included in the analysis. With these parameters we 
can further calculate the column densities for different species 
iV a , and the mean kinetic temperature Tkm- Given the absolute 
intensity of the UV background, the mean gas number density, 
no, and the line-of-sight thickness, L, of the absorber can be 
evaluated as well. 

An advanced version of the MCI completed with the pro- 
cedure of the spectral shape recovering (referred to as MCISS 
hereafter) is iterative and works as follows. At first, some basic 
ionizing spectrum (e.g. HM for corresponding redshift) is taken 
as an initial approximation. To calculate the ionization curves 
needed to carry out the MCI calculations this basic UVB is in- 
serted into the photoionization code CLOUDY (Ferland 1997) 
together with some guesses about relative element abundances 
and metallicity. Then the MCI run is performed with the ioniza- 
tion curves taken as external parameters. If the resulting (recov- 




Fig. 1. Schematic picture of metagalactic ionizing spectrum (solid 
thick line) as predicted by model calculations (Haardt & Madau 1996; 
Fardal et al. 1998). The spectrum is normalized so that J v (hv = 1 
Ryd) = 1. The emission bump at 3 Ryd is caused by reemission of 
Hell Lya, Hell two-photon continuum emission and Hell Balmer 
continuum emission from intergalactic clouds. The positions of ion- 
ization thresholds of different ions are indicated by tick marks. The 
definition of factors {/;} is given in Sect. 2. Dashed line shows the 
spectral shape recovered from metal absorbers at z ~ 3. The contin- 
uum depression between 3 and 4 Ryd (J c — Jgp) is due to He II Lya 
absorption in the IGM (the He II Gunn-Peterson effect). 



ered) relative abundances and the gas metallicity differ from the 
initial guess, the ionization curves are re-calculated with these 
current quantities and new MCI runs are performed till the re- 
covered metallicity and element ratios correspond to those used 
to calculate the ionization curves. If this solution reproduces 
all line profiles within their uncertainties and no other pecu- 
liarities arise (like relative overabundances of elements which 
are well beyond observational and theoretical constraints), then 
we can consider the assumed ionizing spectrum as consistent 
with observations. If, however, the line intensities come over- 
/underestimated and/or relative element ratios are odd, then the 
spectral shape of the ionizing radiation is to be adjusted. 

The spectral shape adjustment procedure itself is developed 
on base of the response surface methodology used in the the- 
ory of experimental design. It includes (i) the parameteriza- 
tion of the spectral shape by means of a set of variables (called 
'factors'), (ii) the choice of a quantitative measure (called 're- 
sponse') to evaluate goodness of a trial spectral shape, and (m) 
the estimation of a direction in the factor space which leads to 
a spectrum with the better goodness. Moving along this direc- 
tion, we come from the initial spectral shape to the one with 
better characteristics concerning the fitting of the observed line 
intensities. Now the MCI calculations are carried out with this 
newly obtained ionizing spectrum and the whole procedure is 
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Fig. 3. Hydrogen and metal absorption lines associated with the z a ^ s = 2.735 system towards HE 2347-4342 (histograms). Synthetic profiles 
corresponding to the recovered ionizing spectrum (dashed line in Fig. 4) are plotted by the smooth curves. The physical parameters are listed 
in Table 1, Col. 2. Bold horizontal lines mark pixels included in the objective function of the MCI procedure. The zero radial velocity is fixed 
at 2 = 2.73576. The central positions of blends are indicated by tick marks. 



repeated till the optimal spectral shape is found, i.e. one which 
allows to reproduce the observed intensities of all lines without 
any physical inconsistencies. 

How well the spectral shape can be recovered depends on 
the number of metal lines involved in the analysis: the more 
lines of different ionic transitions of different elements are de- 
tected in the absorption system, the better constrained is the 
allowable set of the spectral shapes. Any available a priori in- 
formation concerning, e.g., the expected element ratios should 
be also taken into consideration in order to distinguish between 
possible solutions. 

Fig. 1 schematically shows the spectrum of the meta- 
galactic radiation as obtained in model calculations (Haardt & 
Madau 1996; Fardal et al. 1998) and illustrates its parameteri- 
zation by a set of factors. These are: fx - the slope (power law 
exponent) between 1 Ryd and the He II ionization break; /2 - 
the energy at which the He II break starts (point A); fs - the 
depth of the He II ionization break, log( Jb/ </a); fi - the slope 
between A and B; / 5 - the slope between B and C; fa, fa - 
the height and width of the He II re-emission bump. The en- 
ergy of the far UV cut off (point C) when taken above 100 Ryd 
does not effect the fractional ionizations of ions we are interest- 
ing in. In all computations described in the subsequent sections 
this energy and the slope after it are kept fixed at 128 Ryd and 
— 1.5, respectively. 

Studying the UVB at z ~ 3, we found that metagalactic 
ionizing spectra reveal a sharp depression in the intensity be- 
tween 3 and 4 Ryd which was interpreted as manifestation of 



the He II Gunn-Petersen (GP) effect 1 (Paper I). Thus, additional 
factor /g was introduced to describe the depth of the GP depres- 
sion approximated by a straight step. Knowing this depth, Iqp, 
and the continuum level I c at E = 3 Ryd (see Fig. 1), the effec- 
tive He II Lya line opacity due to absorption in the intergalactic 
diffuse gas can be estimated as T cff = ln(/ c // G p)- 

Several absorption systems in the present work have neutral 
hydrogen column densities in the range 10 16 cm -2 < A^(H I) < 
10 17 cm~ 2 , and thus are optically thin in the hydrogen Lyman 
continuum, but can be optically thick in the He II Lyman con- 
tinuum. In Paper I we found that such systems can be described 
with an average ionizing spectrum which reflects both the ex- 
ternal (He II GP absorption) and local effects (strong attenua- 
tion at E > 4 Ryd due to He II continuum absorption in the 
cloud itself and an enhanced peak at E = 3 Ryd due to recom- 
bining emission of He II Lya). An example of such average 
spectrum is plotted by the dashed line in Fig. 1 . 

Fig. 2 shows ionization curves of carbon and silicon ions 
calculated for this spectrum in comparison to those calculated 
for the spectrum of HM (z = 3). The intensity depletion in the 
range 3 < E < 4 Ryd strongly affects ionization fractions of 
carbon leading to enhanced ratios of C Il/C IV and C Ill/C IV at 
log U > — 2.5. The dependence of silicon ion fractions on Iqp 
is rather weak, but they are sensitive to photons with energies 
E < 3 Ryd — an enhanced intensity in this range leads to a 
higher Si iv/Si II ratio. 

1 In the present context we consider the GP effect as being caused 
by combined action of line blanketing and He II Lya absorption in 
continuous intergalactic medium (true GP effect). 
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Fig. 2. Ionization fractions of carbon (upper panel) and silicon (lower 
panel) calculated for metallicity Z = O.OIZq. Thick curves corre- 
spond to the UBV shown by the solid line in Fig. 1, thin curves — to 
the UVB shown by dashed line. 



To complete this section, some words should be said about 
the evaluation of -q. In the systems under study, we observe 
lines of neutral hydrogen, but we do not have spectral ranges 
with the corresponding He II lines. Nevertheless, assuming that 
all lines arise in the same gas and reconstructing the velocity 
and density distributions of the absorbing material from avail- 
able hydrogen and metal line profiles, we are able to estimate 
the column density of single ionized helium as well. The ion- 
ization corrections of He II needed for this estimation are calcu- 
lated using CLOUDY with the recovered UVB spectral shape, 
metallicity and element abundance ratios as inputs. 

3. Analysis of individual absorption systems 

From a wealth of metal absorption systems identified in the 
spectra of HE 2347-4342, Q 1157+3143, and HS 1700+6416 
only a few turned out to be suitable for the UVB recovering. 
These metal absorption-line systems are described in detail in 
the subsections below. 

In general, the uncertainties of the fitting parameters Uq, 
iVn, <J V , &y, and Z a listed in Tables 1, 2 and 3 are about 15%- 
20% (for data with S/N > 30). 

All computations were performed with laboratory wave- 
lengths and oscillator strengths taken from Morton (2003) for 
A > 912 A and from Verner et al. (1994) for A < 912 A. Solar 
abundances were taken from Asplund et al. (2004). Note that 
their solar abundance of nitrogen, N/H = 6 x 10~ 5 , is 1.4 times 
(0.15 dex) lower then that previously reported in Holweger 
(2001). 




Fig. 4. The ionizing background spectrum predicted for z — 3 and 
2 = 2 in HM (thick and thin solid lines, respectively) and the spectra 
restored for z = 2.73 (dashed lines) and z = 1.79 (dotted line). The 
normalization is the same as in Fig. 1. 

The intensity J912 needed to convert the obtained Uq values 
into the gas density was taken from HM. 

3. 1. Absorption systems towards HE 2347-4342 

We used for our analysis the VLT/UVES spectrum of HE 
2347^-342. Observations and data processing are described in 
Shull et al. (2004). 

3.1 .1 . Absorption systems at z = 2.735, 2.739 and 
2.741 

These three absorption systems have a joint column density of 
N(H I) = 2.8 x 10 17 cm -2 (measured from Lya discontinuity 
in the QSO spectrum). The system at z a b s = 2.735 (Fig. 3) ex- 
hibits all silicon and carbon transitions from single to triple ion- 
ized species, N III A989 A (partially blended) and O III A832 A 
(albeit from very noisy part of the spectrum). Expected posi- 
tions of Ovi AA1032, 1037 A are blocked with Lya forest 
absorptions. Test calculations with the HM ionizing spectrum 
ended up with synthetic line profiles which significantly over- 
estimated the observed intensity of Si II and underestimated 
intensities of Si IV and C II. Such behavior of synthetic lines 
along with the procedure of spectral shape adjustment were al- 
ready described in detail in Paper I, Sect.3.1. The final spectral 
shape that we obtained is shown by the dashed line in Fig. 4 
with corresponding physical parameters given in Table 1, Col.2 
and synthetic profiles plotted by the smooth lines in Fig. 3. 

The adjusted ionizing spectrum is much softer at E > 4 
Ryd than the initial spectrum of HM and has a strong emis- 
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sion peak at E = 3 Ryd (needed to reproduce the ob- 
served Si Il/Si IV ratio). At the same time, the observed ratio 
of C Il/C IV cannot be reproduced without a significant drop in 
the intensity between 3 and 4 Ryd. The obtained UVB shape 
resembles that of z abs = 2.9171 system towards HE 0940-1050 
(Paper I, Sect.3. 1) and reflects both local (attenuation at E > 4 
Ryd and enhanced He II Lya emission at E = 3 Ryd) and ex- 
ternal (GP depression at 3 < E < 4 Ryd due to intergalactic 
He II Lya absorption) effects. From the Iqp value the effective 
optical depth of the intergalactic He II Lya absorption at z a bs = 
2.735 is estimated as rgg 11 = 1.45 ± 0.2. The calculated He II 
column density is 7.4 x 10 18 cm -2 which gives 77 = 150. 

Two close absorption systems at z a b s = 2.739 and z a b s = 
2.741 have metal lines either heavily blended or weak and noisy 
leading to the consequence that the observed profiles can be re- 
produced within a broad range of spectral shapes. In any case, 
the UVB recovered from the z a bs = 2.735 system is consistent 
with these systems and delivers reasonable physical parame- 
ters. They are listed in Table 1, Cols. 2 and 3, with correspond- 
ing synthetic profiles plotted in Figs. 5 and 6 (smooth lines). 
Note an extremely low metallicity of [Z] ~ —3.5 in the sub- 
LLS at z a b s = 2.741 along with its huge linear size of L > 250 
kpc. In spite of all uncertainties immanent to this absorber, 
these values can be considered as reliable: the C Ill/C IV ratio is 
similar to that found in the z a b s = 2.735 system (it implies the 
same mean ionization parameter and gas density), but there is 
three times more neutral hydrogen and an order of magnitude 
less C III, C IV and Si IV. As far as we know, this is the lowest 
metallicity ever directly measured in a Lya forest absorber. 

An important result is that a significant difference in 77 be- 
tween the z a b s = 2.739 absorber (77 = 50) and two other ab- 
sorbers (77 ~ 150) stems not from variations in the ionizing 
background (all quantities have been obtained with the same 
ionizing spectrum) but entirely from differences in the gas den- 
sity: n H = 6.0 x 10~ 3 cm" 3 at z abs = 2.739 and (7-9) xl0~ 4 
cm~ 3 at z a bs = 2.741 and 2.735 (see also Fig. 19 below). 

It is also worth noting that these three absorbers are located 
in the vicinity of the He II transmission window at 2.716 < z < 
2.726 marked by 'K' in Figs. 1/ and 5 in Z04. 

In all three systems we measured [Si/C] ~ — 0.1. In or- 
der to interpret this value correctly the following considerations 
should be taken into account. The system at z a bs = 2.735 used 
to restore the UVB is optically thin in H I but opaque in He II 
which supposes that radiation transfer effects can be signifi- 
cant. The average spectrum is appropriate to describe the line 
profiles, but it can produce biased abundances. Computational 
experiments with different spectral shapes of the UVB have 
shown that softening at E > 4 Ryd affects more significantly 
the abundance of silicon making the recovered ratio [Si/C] 
lower. In general, the change of [Si/C] does not exceed 0.2 
dex. This allows to set an upper bound to [Si/C] in the these 
systems: [Si/C] < 0.3. This value is fully consistent with 
[Si/C] measured recently in an extremely metal-poor DLA: 
[C] = -2.83 ± 0.05, [Si] = -2.57 ± 0.05 (Erni et al. 2006), 
and with measurements in metal-poor stars: [Si/C] < 0.3 at 
[Mg] < -3 (Cayrel et al. 2004). 

The systems show a strong underabundance of nitrogen 
compared to the solar value: [N/C] < — 0.5 



3.1 .2. Absorption system at z abs = 1 .7963 

This system was firstly identified through its unusually strong 
Mgll lines (Fig. 7). It also contains absorption lines of Ol 
A1302 A, Cn A1334 A (blended), Sin AA1260, 1193, 1190, 
1526, and 1304 A, Fe n AA2600, 2382 A, Al in AA1854, 1862 
A, Sim A1206 A, Civ A1548 A, and Siiv AA1393, 1402 
A. Unfortunately, positions of N V AA1238, 1242 A coincide 
with Lya forest absorptions, and both lines of the O VI doublet 
fall beyond the available wavelength coverage. From hydrogen 
lines, only a noisy Lya is present and its blue wing is blended. 
However, in spite of this uncertainty the available set of metal 
lines and especially the silicon transitions make it possible to 
estimate the shape of the underlying UVB ionizing radiation. 

Firstly, the O I A1302 A line is observed together with C IV. 
O I traces neutral hydrogen and its ion fraction is practically in- 
dependent on the UVB spectral shape at E > 1 Ryd. The pres- 
ence of O I supposes the ionization parameter U < 10~ 3 , oth- 
erwise the fraction of neutral oxygen becomes too small and the 
resulting overabundance of oxygen to silicon or carbon would 
exceed 1 dex compared to solar values. 

Secondly, all low ionization metal lines exhibit simple pro- 
files and have the same FWHM of 16 km s _1 (fe-parameter 
equals 8.7 km s _1 ), which assumes turbulent motion as the 
main source of the line broadening. Even with this underes- 
timated fe-parameter (thermal component in the broadening of 
the hydrogen line is neglected) the neutral hydrogen column 
density of a few 10 16 cm~ 2 leads already to a significant over- 
estimation of the intensity in the unblended red wing. This 
means that the absorption system is optically thin in the hy- 
drogen continuum and the self-shielding as a reason for the 
simultaneous presence of O I and C IV can be excluded. Thus, 
the only remaining option to observe a pronounced C IV ab- 
sorption at U < 10~ 3 is to have an ionizing spectrum which is 
very hard at E > 4 Ryd. 

Actually, such type of spectrum was recovered by us for the 
absorption system at z a b s = 1-78 towards HE 0141-3932 ex- 
hibiting a similar set of ions (Reimers et al. 2005a). So, it was 
naturally to take that spectral shape as an initial guess for the 
system under study. Test runs have shown that the observed in- 
tensities of silicon lines can be matched only with metallicities 
of 2.5-3Z Q , otherwise the Si III line comes strongly overesti- 
mated for all tried spectral shapes. This behavior is explained 
by a strong dependence of the silicon ion fractions on temper- 
ature which at metallicities near and above solar overrides the 
dependence on spectral shape. 

These considerations allow us to construct the UVB spec- 
trum appropriate for the z a b s = 1 -796 absorption system (Fig. 4, 
dotted line). Due to many uncertainties stated above this spec- 
tral shape is not unique. For example, the power index for the 
slope at 1 < E < 4 Ryd can be higher or lower. This will 
mostly affect the ratio C/Si making it above solar in case of a 
steeper slope. From this point of view a hard slope is prefer- 
able. However, the depth of the He II break at 4 Ryd can be 
constrained: intensities at E = 4 Ryd below O.O6J912 lead 
to a strong overabundance of oxygen to silicon and to carbon 
([Q/Si,C] > 0.2), and ought to be excluded. 
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Table 1. Physical parameters of the z abs = 2.735, 2.739, 2.741 and 1.796 metal absorbers towards HE 2347-4342 derived by the MCISS 
procedure with the recovered UV background spectra shown in Fig. 4. Column densities are given in cm~ 2 
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7 k - 2 739 


7 k - 2 741 


7 k - 1 796 

^abs — 1. 1 7\J 


(i) 


(2) 


(3) 


(4) 


(5) 


C/o 


3.3E-2 


4.4E-3 


4.3E-2 


2.0E-3 




2.3E20 


7.3E19 


5.9E20 


2.0E17 


ov, km s 1 


25.0 


24.0 


21.0 


15.0 


°V 


0.6 


0.5 


0.7 


0.6 


Zc 


2.4E-6 


1.2E-6 


7.3E-8 


7.4E-4 


Zn 


<2.0E-7 


< 8.0E-8 






Zo 


5.0E-6 






1.3E-3 


Z-Mg 








(1.4-2.5)E-4 


Za\ 








1.1E-1 


Zsi 


4.3E-7 


1.5E-7 


1.0E-8 


9.5E-5 


Zys 








2.3E-5 


[Zc] 


-2.00 


-2.30 


-3.53 


0.48 


[Zn] 


< -2.5 


< - 2.87 






[Zo] 


-2.00 






0.45 


[Zm s ] 








0.60-0.85 


[Zm] 








0.6 


[Zsi] 


-1.90 c 


-2.34 c 


-3.53 c 


0.46 


[Z Fc ] 








-0.08 


N(H I) 


(4.9±0.5)E16 


(8.6 ± 1.0)E16 


(1.4±0.2)E17 


(9.1±2.0)E15 


JV(Oi) 








(1.0±0.2)E13 
S^ED 6 


N(C II) 


(7.1±0.5)E12 


(4.1±0.4)E12 


(6.2±0.6)E11 


N(Mg II) 








(0.9-1.6)E13 


JV(Sin) 


(1.6 ± 0.2)E12 


(5.5±0.5)E11 


(1.0±0.2)E11 


(1.4±0.1)E13 


7V(Fe II) 








(9.4±0.5)E11 


JV(Cra) 


(3.7±0.4)E14 


6.9E13 6 


3.0E13 b 




N(Nm) 


< 2.5E13 


< 4.6E12 






iV(Om) 


(7.5±2.0)E14 








JV(A1 in) 








(7.4±0.7)E11 


JV(Silll) 


(2.5 ± 0.2)E13 


6.9E12 6 


\m\2 h 


(4.9 ± 1.0)E12 


N(C IV) 


(1.2 ± 0.1)E14 


(4.4 ± 0.5)E12 


(9.5 ± 1.5)E12 


(1.2 ± 0.1)E13 


iV(Siiv) 


(2.0±0.2)E13 


2.0 ± 0.2)E12 


(1.1 ±0.1)E12 


(2.9±0.3)E12 


(T),K 


2.9E4 


2.0E4 


3.2E4 


5.5E3 


riH, cm 


9.0E^1 


6.0E-3 


7.0E-4 


2.2E-2 


L, kpc 


90 


4 


280 


0.004 


7V(He II) 


7.4E18 6 


4.3E18 6 


2.0E19 6 


2.0E16 ft 


V 


150 ± 20 


50 ±5 


140 ± 20 


2 ± 2 


tE§> u (at 3 Ryd) 


1.45 ±0.2 










a Z x = N x /Nn; [Z X ] = \og(N x /NH) - log{N x /N H ) @ 

^calculated using the velocity-density distributions estimated from hydrogen and metal profiles 
^actual abundance may be 0.2 dex higher — see the last paragraph in Sect. 3.1.1 



Physical parameters estimated with the UVB from Fig. 4 
are given in Table 1, Col. 5. Strong relative underabundance of 
iron (three times solar values for carbon, oxygen and silicon 
and only 0.8 for iron) can be considered as reliable since it was 
reproduced for all tried spectral shapes. Aluminium is slightly 
overabundant compared to C, O and Si (4 times solar), but the 
accuracy of this estimation is not high enough due to uncertain- 
ties in the recombination coefficients used in CLOUDY. As for 
magnesium, a spread of its abundance (4-7 solar) is explained 
by the saturation of Mg II A2796 A and possible contamination 
of Mg II A2803 A with a telluric line. 

3.2. Absorption systems towards Q 1 157+3143 

Absorption spectrum of this quasar was obtained with 
Keck/HIRES. Details on the data reduction are given in 



Kirkman & Tytler (1999). Additionally, this bright quasar was 
observed with HST/STIS to study He II Lya absorption in the 
intergalactic matter (Reimers et al. 2005b). 

3.2.1 . Absorption system at z abs = 2.944 

This absorption system (Fig. 8) lies very close to the quasar 
(z om = 2.97, Lya forest starts at z = 2.999) and could be as- 
sociated. However, Q 1 157+3143 does not show proximity ef- 
fect. On the contrary, the line density in its vicinity is unusually 
high (Ganguly et al. 2001). A set of ions detected in the z abs = 
2.944 system is also typical for an intergalactic absorber: asso- 
ciated systems are usually highly ionized and show pronounced 
lines of O VI, N V, and C IV along with weak C III, Si III, and 
S IV, whereas the system under study reveals C III and Si III 
lines stronger than C IV, Si IV, and no N V absorption. At the 
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Fig. 6. Same as Fig. 3 but for the z a h s = 2.741 system (HE 2347^-342). The zero radial velocity is fixed at z = 2.74097. The estimated physical 
parameters are listed in Table 1, Col. 4. 



expected position of the Ovi A1031 A line only a weak ab- 
sorption is seen (its identification as O VI is ambiguous since 
Ovi A1037 A is blended). 



From the H I continuum absorption in the observed spec- 
trum the neutral hydrogen column density of this system (to- 
gether with its neighbor at z a h s = 2.939) is estimated as (1.5 ± 
0.2) x 10 17 cm~ 2 . Such a high hydrogen content together with 
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Fig. 7. Same as Fig. 3 but for the z a hs = 1.796 system (HE 2347-4342). Synthetic profiles (smooth curves) were calculated with the recovered 
ionizing spectrum shown by the dotted line in Fig. 4. The zero radial velocity is fixed at z = 1.7963. The estimated physical parameters are 
listed in Table 1, Col. 5. Synthetic profiles of Mg II are calculated for the abundance 2.5 x 1CT 4 (see text). Blends are indicated by tick marks. 



rather weak intensities of the metal lines point to a very low 
metallicity of the absorbing gas which is another argument 
against the associated nature of the z a b s = 2.944 system. 

Test calculations with the HM ionizing spectrum confirmed 
the metallicity of 0.001.Z©. However, with this spectrum the 
calculated intensity of the C III line in the range — 70 < v < 70 
km s _1 was significantly less than observed. Since the position 
of this line is in the Lya forest, possible blends should be taken 
into account. We do not find any metal absorption contaminat- 
ing the C III line. The only candidate for blending could be a 
hydrogen absorption at z a b s = 2.17. If we assume that observed 
profile is mainly caused by this hydrogen absorption, then the 
kinetic temperature of the gas estimated from the line width is 
Tkin < 32000 K. This implies that the ionization parameter U 
is less than 0.02 (for metallicity \Z\ < —1.5, since there are 



no metal lines detected at z a b s = 2.17) and the linear size of 
the putative H I absorber does not exceed 0.5 kpc. Such small 
cloud would hardly survive in the IGM. Another argument in 
favor of C III is the comparison with the Si III line: both lines 
have similar widths (FWHM = 40 km s _1 and 36.5 km s _1 ) 
and the same asymmetry. Thus, we can assume with high prob- 
ability that the absorption at —70 < v < 70 km s _1 is indeed 
caused by C III from the z a b s = 2.944 system. In this case only 
a spectrum with a significant flux depression between 3 and 4 
Ryd makes it possible to fit adequately all observed lines. 

The reconstructed spectral shape of the UVB is shown in 
Fig. 9 by the dotted line, the corresponding physical parame- 
ters are given in Table 2, Col. 2, and the synthetic profiles are 
plotted by the solid lines in Fig. 8. 
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Fig.8. Same as Fig. 3 but for the z a b s = 2.944 system towards Q 1157+3143. Synthetic profiles (smooth curves) were calculated with the 
recovered ionizing spectrum shown by the dotted line in Fig. 9. The zero radial velocity is fixed at z — 2.9444. The estimated physical 
parameters are listed in Table 2, Col. 2. Blends are indicated by tick marks. 



Here we find again a spectrum similar to that from the 
z abs = 2.9171 system (HE 0940-1050, Paper I): its shape re- 
flects both the external (GP depression) and local effects (an 
enhanced emission peak at 3 Ryd and a sharp break at 4 Ryd) 
arising from the reprocession of the incident ionizing radiation 
by the absorbing cloud itself. 

The Hell Lya opacity is estimated as Tqp 11 = 2.5 ± 0.2. 
The predicted value for the He II column density is 8.5 x 10 18 
cm~ 2 which gives r] = 85. 

The ratio [Si/C] = 0.06, but as mentioned in Sect. 3.1.1 the 
actual value may be 0.2 dex higher: [Si/C] < 0.26. 

A significant relative underabundance of nitrogen, [N/C] < 
—0.3, is also worth noting for this system. 

3.2.2. Absorption system at z abs = 2.939 

This system presents metal lines of both low (C II, Si II) 
and high (Civ, Si IV) ionization stages (Fig. 10). Lines Ovi 
AA1031, 1037 A are blended with strong Lya forest absorp- 
tions. The spectrum of HM turned out to be inconsistent with 
the observed line intensities, the adjusted spectral shape is 



shown in Fig. 9 by the short-dashed line. The corresponding 
physical parameters are given in Table 2, Col. 3, and the syn- 
thetic profiles are plotted by the solid lines in Fig. 10. 

Although shifted by only Az = 0.0045 (Av = 344 
km s _1 ) from the z a bs= 2.944 system, the present absorber 
reveals significantly different physical conditions: metallicity 
two order of magnitude higher along with a lower ionization 
parameter (higher gas density) and a very small (for an inter- 
galactic absorber) linear size below 1 kpc. 

The z a b s = 2.939 system has optical depth ~ 1 in the He II 
Lyman continuum. Thus, the attenuation of the incident radia- 
tion is probably not very pronounced and the recovered spec- 
tral shape can represent in this case the external UVB. It is seen 
from Fig. 9 that spectrum is softer at E > 4 Ryd than the model 
spectrum of HM and has a more pronounced He II Lya emis- 
sion peak at E = 3 Ryd. This is probably a consequence of a 
high number density of clouds registered at z < 3 along the 
sight line towards Q 1 157+3143. 

The predicted value of rj = 45 ± 5 is two times lower than 
r] in the z abs = 2.944 system and reflects both the higher gas 
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Table 2. Physical parameters of the z a bs= 2.944, 2.939, and 2.875 metal absorbers towards Q 1157+3143 derived by the MCISS procedure 
with the recovered UV background spectra shown in Fig. 9. Column densities are given in cm~ 2 



Pa ra TTi ptpr a 

L CIL ell 111^ L\^L 


7 k - 2 944 


7 u - 2 939 


7 k - 2 875 


(1) 


(2) 


(3) 


(4) 


U 


1.6E-2 


8.8E-3 


6.9E-2 


N H 


1.4E20 


8.3E18 


2.4E20 


er v , km s 1 


20.0 


22.5 


25.0 


°y 


0.6 


0.6 


0.5 


Zc 


3.0E-7 


2.3E-5 


1.1E-6 




<3.0E-8 


<1.0E-6 




Z S i 


4.7E-8 


3.7E-6 


1.6E-7 


[Zc] 


-2.90 


-1.02 


-2.36 


[^n] 


< -3.3 


< -1.8 




[*Si] 


-2.84 d 


-0.94 


-2.34 d 


iV(Hi) 


(1.1 ±0.2)E17 


1.2 ±0.1)E16 


(1.6 ± 0.4)E16 


N(C II) 


(1.4±0.2)E12 


(9.0±0.5)E12 


(2.0 ± 0.3)E12 


/Y(Sin) 


(2.5±0.2)E11 


(1.5±0.2)E12 


(1.8±0.2)E11 


N(C m) 


(3.7±0.6)E13 i ' 


(1.8±0.2)E14 


(1.9 ± 0.2)E14 


N(N m) 


< 5.8E12 


< 7.5E12 




/V(Siin) 


(3.4 ± 0.4)E12 b 


2.0E13 C 


4.2E12 C 


7Y(C IV) 


(3.2±0.8)E12 


(1.1 ±0.3)E13 


Mn + n d'iFn 


TV(Siiv) 


(1.7±0.3)E12 


(8.0±0.4)E12 


(6.8±0.5)E12 


(T>,K 


2.7E4 


1.8E4 


3.7E4 


n-H, cm ' 


2.0E-3 


3.0E-3 


4.0E^1 


£, kpc 


27 


0.86 


200 


iV(He 11) 


8.5E18 C 


5.4E17 C 


7.0E18 C 




85 ± 10 


45 ±5 


440 ± 40 


r£| n (at3Ryd) 


2.5 ±0.3 


2.5-3 



a Z K = N K /N H ; [Z X ] = log(N x /N H ) - \og(N K /N H ) e ; partially blended; 

c calculated using the velocity-density distribution estimated from hydrogen and metal profiles 

d actual value may be 0.2 dex higher — see the last paragraph in Sect. 3.1.1 



3.2.3. Absorption system at z abs = 2.875 

Hydrogen and metal lines identified in this system are shown in 
Fig. 11. The C iv A1548 A line falls in a gap between echelle 
orders, Nlll A989 A and Ovi AA1031,1037 A are severely 
blended. Note an extremely weak Si II A1260 A line which is 
nevertheless very well detected owing to high S/N = 150 in this 
wavelength range. 

Test calculations with the HM ionizing spectrum showed 
carbon content [C/H] = -2.3 associated with a very high rel- 
ative overabundance of silicon, [Si/C] ~ 0.8. The calculated 
intensity at the expected position of C III came significantly un- 
derestimated, but blending cannot be excluded (the only candi- 
date — H I Lya at z a b s = 2. 1 15). 

The hydrogen column density is not well defined because 
of high noise at the locations of Ly-7 A926 A and Ly-8 A923 
A, but, in any case, it does not exceed 3 x 10 16 cm~ 2 . 

The obtained overabundance of silicon points to the ioniz- 
ing spectrum much softer at E > 4 Ryd than the HM spectrum 
that we used. A chance that high [Si/C] arises from contribu- 
tion of the hypothetical Pop III stars with masses M > 100M Q 
(pair-instability SNe) seems quite unprobable since such over- 
abundance has never been observed even in the extremely 
metal-poor clouds with Z < 0.001Z o (see Sect. 3.1.1). A soft 
spectrum can be produced either by the He II continuous ab- 
sorption (both in the IGM and inside the cloud itself) or by the 
input of radiation from soft sources like stars. In the present 
case the first option is preferable taking into account the ob- 
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Fig. 9. Same as Fig. 4 but for the systems at z abs = 2.944, 2.939, and 
2.875 towards Q 1157+3143. 



density and the ionizing background which is harder at E > 4 
Ryd (see Fig. 9). 

Noteworthy is again a very strong relative underabundance 
of nitrogen, [N/C] < -0.8. 
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Fig. 10. Same as Fig. 3 but for the z abs = 2.939 system towards Q 1157+3143. Synthetic profiles (smooth curves) were calculated with the 
recovered ionizing spectrum shown by the short-dashed line in Fig. 9. The zero radial velocity is fixed at z = 2.93987. The estimated physical 
parameters are listed in Table 2, Col. 3. Blends are indicated by tick marks. 




Fig. 11. Same as Fig. 3 but for the z abs = 2.875 system towards Q 1157+3143. Synthetic profiles (smooth curves) were calculated with the 
recovered ionizing spectrum shown by the long-dashed line in Fig. 9. The zero radial velocity is fixed at z — 2.87567. The estimated physical 
parameters are listed in Table 2, Col. 4. Blends are indicated by tick marks. 



servational evidence of large amount of He II at z a b s — 2.87 
towards Q 1 157+3143 (Reimers et al. 2005b), and a high value 
of the ionization parameter U (log(Si iv/Si II) ~ 1.6) which 
supposes 7] > 150 (see Fig. 19) which assusmes, in turn, that 
the cloud is opaque in He II. 



Unfortunately, high noise in the C II A1334 A profile and 
blending of Si III hamper the accurate reconstruction of the 
ionizing background. Calculations with the spectrum from the 
z a bs = 2.944 system lead to the silicon overabundance [Si/C] ~ 
0.2-0.3. Taking into account that the average spectrum gives 
[Si/C] by < 0.2 dex lower than the real value, this translates 
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Fig. 12. Same as Fig. 3 but for the 2 a b s = 2.568 system towards HS 1700+6416. Synthetic profiles (smooth curves) were calculated with the 
recovered ionizing spectrum shown by the dotted line in Fig. 13. The zero radial velocity is fixed at z = 2.56826. The estimated physical 
parameters are listed in Table 3, Col. 2. 



to [Si/C] « 0.4 — 0.5 which is slightly higher than measured 
elsewhere (Sect. 3.1.1). A UVB with a lowered by 0.2 dex de- 
pression between 3 and 4 Ryd (Fig. 9, long-dashed line) gives 
[Si/C] ~ (and reproduces the observed line intensity at the 
position of C III). These two shapes constrain the acceptable 
range. The physical parameters obtained with the latter spec- 
trum are listed in Table 2, Col. 4. The corresponding synthetic 
profiles are plotted by the solid lines in Fig. 11. The calculated 
column density of He II, A^(He II) = 7.0 x 10 18 cm~ 2 , sup- 
poses that the cloud is optically thick in the He II continuum 
and, thus, the restored spectral shape is again a superposition 
of the local and external effects. 

The mean optical depth of Hell Lya, Tq£ h = 2.5-3, is 
comparable to r^p 11 = 2.09 ± 0.1 estimated at z a b s = 2.87 by 
Reimers et al. (2005b). 

We note the extremely high value of r\ — 440 ± 40 which 
is reliably reproduced in all test runs in spite of uncertainties in 
both the neutral hydrogen content and the spectral shape of the 
ionizing radiation. This value can be explained by the softness 
of the effective ionizing radiation field and by the low gas den- 
sity in the absorbing cloud (cf. with r\ = 150 at z a b s = 2.944 
produced by the ionizing spectrum of similar softness). 

3.3. Absorption systems towards HS 1700+6416 

This one of the brightest known high redshift quasar was ob- 
served many times with many instruments, both ground- and 
satellite-based. Reimers et al. (1992) and Vogel & Reimers 
(1995) firstly reported on 7 strong metal absorption systems 
identified in its low-resolution combined optical and UV spec- 
tra. Later on, Davidsen et al. (1996) measured opacity of the in- 
tergalactic He II Lya using far UV spectrum obtained with the 
HUT. Simcoe et al. (2002, 2006) studied physical conditions in 
the metal absorption systems using high-resolution Keck tele- 
scope data. A comprehensive analysis of both H I and He II 
(FUSE) spectra was performed by F06 and Fechner & Reimers 
(2006). For our reconstruction of the underlying UVB spectra 



we selected systems with many different ionic lines, that are in 
the ionization equilibrium and do not show self-shielding ef- 
fects 2 . Observational data stem from the Keck/HIRES (AA = 
3150-9000A)andHST/STIS (AA = 1155-1700 A) spectra. 
Details on data reduction are given in Simcoe et al. (2002) and 
Fechner et al. (2006b). 

3.3.1 . Absorption system at z ahs = 2.568 

Hydrogen and metal lines identified in this system are shown 
in Fig. 12. Test calculations with the ionizing spectrum of HM 
revealed that the red wing of H I Lya is inconsistent with the 
assumption of constant metallicity inside the absorber. In the 
following, the red wings of hydrogen lines are reconstructed 
from the density-velocity distributions derived from metal lines 
under the assumption of constant metallicity. 

Calculations with the HM spectrum produced also strongly 
underestimated intensity of C III. Blending of C III seems rather 
unprobable since the line is weak, narrow and its asymmetric 
profile is similar to those of C IV and O VI. Thus the shape of 
the HM ionizing spectrum was modified. The result is shown 
in Fig. 13 by the dotted line. The corresponding physical pa- 
rameters are given in Table 3, Col. 2, and the synthetic profiles 
are plotted by the solid lines in Fig. 12. 

The mean optical depth of the intergalactic He II Lya ab- 
sorption decreases at redshift z a b s = 2.56 to the value r^js 11 = 



2 Since metal systems towards HS 1700+6416 are well known, the 
reader may be interested in more details why other systems were not 
included in the UVB recovering procedure: z = 2.578 — two over- 
lapping systems with different metallicities which cannot be clearly 
separated; z = 2.315 — probably a non-equilibrial system with overi- 
onized O VI, may be associated with an AGN/QSO close to the line of 
sight; z = 2.168 — the ratios Si IV/Si II ~ 40 and Si IV/C IV ~ 0.16 
point out to a very soft UVB (cf. data for the z = 2.875 system in 
Table 2), however, uncertainties in C II and Si III (blending ?) and the 
absence of C III A977 A do not allow us to reconstruct the actual spec- 
tral shape; z — 1.724 — most metal lines are contaminated by the 
Lya forest absorptions. 
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Table 3. Physical parameters of the z abs = 2.568, 2.438, 2.433, 2.379, and 1.845 metal absorbers towards HS 1700+6416 derived by the MCISS 
procedure with the recovered UV background spectra shown in Fig. 13. Column densities are given in cm~ 2 
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1.2 ± 0.2. Note again a significant underabundance of nitrogen, 
[N/C] < -0.3. 



3.3.2. Absorption systems at z ahs = 2.438 and z abs = 
2.433 

Both absorption systems (Figs. 14 and 15) are parts of a large 
absorption complex stretching over 1000 km s _1 (see Fig. 6 in 
Simcoe et el. 2006) and represent dense clumps embedded in 
a highly ionized gas traced by O VI. While this highly ionized 
gas could be still non-equilibrial, the systems under considera- 
tion have already reached the ionization equilibrium with the 
background radiation owing to their higher gas density and, 
hence, a shorter recombination time. Note that the whole ab- 
sorbing complex resembles the z a b s = 1.8073 complex towards 
HS 0747+4259 described in Sect. 4.1 in Paper II. 



From low-resolution spectrum, the neutral hydrogen col- 
umn density of the whole complex is estimated as 4.0 x 10 16 
cm~ 2 (Vogel & Reimers, 1995). Test calculations with the HM 
ionizing spectrum resulted in an underestimated intensity of 
C II A1334 A along with a strong overestimation of Si II A1193 
A. The adjusted spectral shape is shown in Fig. 13 by the short- 
dashed line. The calculated physical parameters are listed in 
Table 3, Cols. 3 and 4, and the corresponding synthetic profiles 
are shown by the solid lines in Figs. 14 and 15. 

Although the spectral shape differs from that recovered for 
the z a bs = 2.568 system, the mean He II optical depth is the 
same: r^p 11 = 1.2 ± 0.2. This value is to be compared with 
r Gr> n = 1-0 ± 0.07 measured at z — 2.4 in the low resolution 
far UV spectrum of HS 1700+6416 by Davidsen et al. (1996), 
and with rfg 11 = 0.74 ± 0.34 at z = 2.45 measured from 
FUSE data by F06. It should be noted, however, that latter value 
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Fig. 14. Same as Fig. 3 but for the 2 a b s = 2.438 system towards HS 1700+6416. Synthetic profiles (smooth curves) were calculated with the 
recovered ionizing spectrum shown by the short-dashed line in Fig. 13. The zero radial velocity is fixed at z — 2.4386. The estimated physical 
parameters are listed in Table 3, Col. 3. Blends are indicated by tick marks. 
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Fig. 15. Same as Fig. 3 but for the z a b s = 2.433 system towards HS 1700+6416. Synthetic profiles (smooth curves) were calculated with the 
recovered ionizing spectrum shown by the short-dashed line in Fig. 13. The zero radial velocity is fixed at z — 2.4330. The estimated physical 
parameters are listed in Table 3, Col. 4. 



can be underestimated due to problems with zero level flux in 
the FUSE spectrum (see Sect. 2 in F06). 

Another point worth noting is the variation we find in 77: 
77 = 45 at z abs = 2.433 and 77 = 100 at z abs = 2.438 (Table 3), 
which occurs because of the density differences in the clouds 
under study. 

The value of r/ at z a b s = 2.438 is close to the 77 estimated 
for the previous z a b s = 2.568 system although the gas density 
in these two absorbers is quite different. This result is a conse- 
quence of softer UVB at z abs = 2.438 (see Fig. 19 in Sect. 4). 

3.3.3. Absorption system at z ahs = 2.379 

A set of metal lines observed in this system (Fig. 16) is identical 
to that from the z a b s = 2.568 absorber described in Sect. 3.3.1. 
Unfortunately, a noisy profile of C III hampers the accurate re- 
construction of the spectral shape of the underlying radiation. 
Although this line lies deep in the Lya forest, it is very narrow 



(FWHM — 22.5 km s -1 ) and cannot be a hydrogen line. On 
the other hand, no metal candidate for possible blending was 
found. Thus, the observed absorption is with high probability 
due to C III. However, the spectrum of HM underestimates the 
observed intensity of C III far above the noise level, whereas the 
spectrum derived for the z a b s = 2.433 system underestimates 
Si III. At the same time, the spectral shape from z a b s = 2.568 
system turned out to be quite consistent with all observed in- 
tensities. The physical parameters obtained with this UVB are 
given in Table 3, Col. 5 with the corresponding synthetic pro- 
files shown by the solid lines in Fig. 16. 

3.3.4. Absorption system at z abs = 1 .845 

This system reveals one hydrogen line (Lya) and a lot of met- 
als in different ionic transitions (Fig. 17) which allow us to 
probe the underlying ionizing background. The column density 
of neutral hydrogen was estimated from low resolution spec- 
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Fig. 16. Same as Fig. 3 but for the 2 a b s = 2.379 system towards HS 1700+6416. Synthetic profiles (smooth curves) were calculated with the 
recovered ionizing spectrum shown by the dotted line in Fig. 13. The zero radial velocity is fixed at z — 2.3799. The estimated physical 
parameters are listed in Table 3, Col. 5. 



cm civ 

Kill ' NV 

OIII OIV ' OVI 




1 10 E, Ryd 

Fig. 13. Same as Fig. 4 but for the systems at z a b s = 2.568, 2.433, 
2.438, 2.379, and 1.845 towards HS 1700+6416. 



trum by Vogel & Reimers (1995) as 5.6xl0 16 cm -2 , i.e. this 
system is optically thin in the hydrogen continuum. It is thin in 
the He II Lyman continuum as well: the attenuation of the inci- 
dent radiation due to continuous He II absorption in the cloud 
itself leads to an increasing Si IV/C IV ratio, but we find in this 
system Si IV/C IV = 0.088 which is well below the values mea- 
sured in the systems with local attenuations at z a b s = 2.735, 
2.944, and 2.875 (Tables 1 and 2). Assuming optical depth in 
He II continuum below 1 one obtains from A^He II) = 77- A^(H I) 
< 5 x 10 17 cm -2 the vlue of 77 ~ 10 which assumes a low ion- 
ization parameter (high gas density) and a hard ionizing spec- 
trum. 

A hard ionizing continuum is also required by the presence 
of N V, otherwise we obtain [N/C] > 0. Such overabundance of 



nitrogen is predicted at high over solar metallicities by chem- 
ical evolution models (e.g. Hamann & Ferland 1999), but at 
Z ~ Z Q still deficit of nitrogen, [N/C] < 0, is measured (e.g. 
Srianand & Petitjean 2001; Reimers et al. 2005a; Jenkins et 
al. 2005). In principle, we cannot exclude a possibility that the 
present system is out of equilibrium and N V is still overionized 
(i.e. not in the ionization equilibrium with the background ra- 
diation). The O VI lines could clear up the actual physical con- 
ditions but they fall beyond the observable range. However, we 
do not observe any velocity shifts between the N V lines and all 
other ionic transitions which could indicate a multiphase na- 
ture of the absorbing gas. This makes 'normal' (equilibrium) 
conditions preferable. 

Another very important clue to the physical properties 
of the absorbing gas deliver silicon ions: the observed ra- 
tios of Si iv/Si III ~ 1 and Siiv/Sill ~ 10 can be realized 
only at metallicities near solar values. As already explained in 
Sect. 3.1.2, this stems from a strong dependence of silicon ion 
fractions on the temperature which at high metallicities over- 
rides dependence on the spectral shape of the UVB. 

Taking together, these considerations allow us to design a 
possible shape of the ionizing continuum. It is shown in Fig. 13 
by the long-dashed line with the derived physical parameters 
given in Table 3, Col. 6 and synthetic profiles plotted by the 
solid lines in Fig. 17. The red wing of a sole hydrogen line is 
inconsistent with the assumption of constant metallicity inside 
the absorber and we calculated it using the density-velocity 
distribution obtained from metal lines. As mentioned above, 
the available silicon lines require metallicity close to solar. 
However, with Z ~ 0.5Z Q we obtained 7Y(Hl) = 2.5 x 10 16 
cm~ 2 which is already half that estimated from the Lyman 
edge in low resolution spectrum. A possible explanation of 
this discrepancy may be that we observe a high metallicity 
(Z ~ Z Q ) blob with 7Y(Hl) - (1-2) xlO 16 cm" 2 embedded 
in a cloud or a chain of clouds with a joint column den- 
sity N(H I) ~ (3-4) xlO 16 cm~ 2 . To confirm this assumption, 
higher order H I Lyman series lines are needed. 
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Fig. 17. Same as Fig. 3 but for the 2 a b s = 1.845 system towards HS 1700+6416. Synthetic profiles (smooth curves) were calculated with the 
recovered ionizing spectrum shown by the long-dashed line in Fig. 13. The zero radial velocity is fixed at z = 1.8450. The estimated physical 
parameters are listed in Table 3, Col. 6. Blends are indicated by tick marks. 



Because of this uncertainty of iV(H I) in the metal-bearing 
cloud and unknown element ratios, the spectral shape plotted 
in Fig. 13 is not unique. However, it delivers quite reasonable 
and self-consistent physical parameters and, thus, can be con- 
sidered as a plausible solution. In any case, the spectral indices 
a > —1.2 in the range 1 < E < 4 Ryd are favored since 
for lower indices one would obtain C/Si above solar which has 
never been observed and is not predicted theoretically. 



Taking together, the high metallicity of the blob and hard 
ionizing spectrum allow to suggest the presence of a QSO/AGN 
close to the line of sight. 



4. Discussion 

The spectral shape of the background ionizing radiation in the 
energy range 1 < E < 10 Ryd is reconstructed from optically- 
thin metal absorption-line systems observed in QSO spectra at 
z « 1.80 and 2.38 < z < 2.94. The spectral shape of the UVB 
does fluctuate, however these fluctuations are ruled by different 
sources at two considered redshift intervals. 

In the range 2.38 < z < 2.94, the spectral shape is mostly 
affected by radiative transfer processes in the clumpy IGM in- 
volving absorption in both the He II Ly-a and Lyman contin- 
uum. We confirm our result from Paper I: the UVB at 2.4 < 
z a bs < 3.0 demonstrates the intensity depletion between 3 and 
4 Ryd caused by the absorption of the Hell Lya photons in 
the IGM (He II Gunn-Petersen effect). None of the recovered 
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UVB shapes shows features which can be attributed to the con- 
tribution of radiation from soft sources like starburst galaxies. 
Such contribution would lead to a very steep slope at E > 1 
Ryd (a < -2.0 at 1 < E < 3 Ryd) and reduced intensity 
at 3 Ryd (see Figs. 1 and 3 in Giroux & Shull, 1997) whereas 
the metal-line systems considered here require relatively hard 
(a > —1.6) spectral shapes with pronounced peaks at 3 Ryd 
combined with strong attenuation at E > 4 Ryd — features un- 
ambiguously related to the absorption in the He II Lyman con- 
tinuum and subsequent re-emission of recombinant He II Ly-a. 

From the depth of the GP trough in the recovered UVB 
spectra we can estimate the effective He II Lya opacity of the 
IGM, Tqp 11 . Fig. 18 shows the obtained redshift evolution of 
r Gp n (filled symbols) together with the available directly mea- 
sured He II Lya opacities (open symbols) along different sight- 
lines. The concordance between the recovered values and the 
direct measurements is good although our estimations are sys- 
tematically higher by 20-40% (~ l-2er). This difference can 
be explained by the following reasons. Firstly, the approxima- 
tion of the continuum depression between 3 and 4 Ryd by a 
straight step does not take into account a winding shape of the 
actual GP trough. Secondly, the depth of this step is determined 
mostly from the analysis of carbon lines C II, C III, and C IV 
(see Sect. 2). However, the ionization potential of Cm (3.52 
Ryd) is close to the onset of the Hell Ly/3 forest (3.56 Ryd) 
and an additional GP absorption in He II Ly/3 can affect the 
fraction of this ion and the subsequent C IV. These may lead 
to the overestimation of r^p 11 by 0.1-0.2 in our method. On 
the other hand, estimations of r^p 11 from FUSE data may give 
too small values due to problems in background subtraction at 
short wavelengths corresponding to z < 2.55 (see Sect. 2 in 
F06). 

Analysis of the close absorption systems (z a b s = 2.735 and 
2.739 towards HE 2347^342 and z abs = 2.433 and 2.438 to- 
wards HS 1700+6416) revealed two-three times fluctuations in 
r\ between them. These systems share the same UVB, and the 
fluctuations in r\ are caused entirely by the variations in the ab- 
sorber gas density. 

Now compare the z a b s = 2.433 and 1.845 systems. Both of 
them show similar U and rin, but their r\ values differ by about 
3 times (Table 3). In this case we observe the effect of different 
spectral softness (see Fig. 13). 

Fig. 19 shows the dependence of r? mod e] = ZhoYhcIi/Yhi 
on the ionization parameter U calculated with CLOUDY for 
optically thin case assuming different ionizing spectra. In case 
of a constant gas density inside the absorber, r7 mo d c i = ??■ 
When the gas density fluctuates within the absorber, rj is de- 
fined as a density weighted mean of the rj{U) values, thus be- 
ing lower than r? mo dci corresponding to the mean ionization 
parameter Uq (see eqs. [24,25] in LAK). It is seen that in the 
range —3 < log(£7) < — 1 relevant for optically thin high red- 
shift intergalactic absorbers the 77 mo dercurve demonstrates a 
large gradient which becomes steeper with the softening of the 
UVB at E > 4 Ryd. Moreover, the interval of ?y m odci values 
corresponding to —3 < \og(U) < — 1 also expands. Thus, in 
general, the scatter in the measured values of r\ can arise from 
the combined action of differences in the mean gas density be- 
tween the absorbers and variable softness of the ionizing radi- 



ation due to absorption in the He II Lyman continuum. We can 
also predict that the scatter in r\ should decrease with decreas- 
ing z due to progressive hardening of the UVB spectrum. 

Fig. 19 explains also the correlation between low Hi col- 
umn densities and high rj values: weak absorbers tend to have 
lower gas densities and thus are shifted to the right part of the 
rj(U) curve. This is also one of the reasons why r/ in voids 
are systematically higher than in the filaments: by definition, 
voids consist of weak H I absorbers. Another reason is a pos- 
sible shielding of the UV radiation in voids by the surround- 
ing clouds which can soften the ionizing spectra. This is il- 
lustrated by the z abs = 2.735, 2.739 and 2.741 absorbers to- 
wards HE 2347^1-342 located at the boundary of the void at 
2.718 < z < 2.727 (Fig. 5 in Z04). These absorbers belong 
probably to the filamentary structure and form a layer with a 
linear size of ~ 300 kpc 3 . Table 2 from Z04 shows that the 
T] values in the void at 2.718 < z < 2.727 are significantly 
higher than in the void at 2.816 < z < 2.823. The latter 
has no strong Hi absorbers at its boundaries and lies close 
to the quasar HE 2347-4342, i.e. the UVB in this void can 
be quite hard. However, over large ranges in z the correlation 
"low N(H I) - high 77" [N(Ui) = 10 12 - 10 13 cm" 2 ] should 
be weak due to superposition of UVBs with variable softness: 
harder UVBs give a significantly lower 77 especially for low 
gas density (and, in general, low column density) absorbers, 
thus effectively lowering the left side of the " N(H 1) - 77" dis- 
tribution. As for stronger H I absorbers [N(H I) = 10 13 — 10 15 
cm~ 2 ], they have usually higher gas densities (lower U values) 
and according to Fig. 19 their 77 values are less sensitive to the 
variations of UVB. Thus, the correlation "high 7V(Hl) - low 
77" should hold over all redshifts. This is in concordance with 
direct measurements of 77 from the H I/He II Lya forest in F06 
(see their Fig. 10) and Fechner & Reimers (2006). 

As for UVBs recovered at z ~ 1.8, they are very hard 
at E > 3 Ryd — much harder than predicted by the HM 
model (Fig. 20). The reason for this discrepancy was already 
given in our Paper II: HM used a biased statistics of the Lya 
clouds at z a b s < 1.8 with too many strong Lya absorbers 
leading to a significant underestimation of the UVB inten- 
sity at E > 3 Ryd. From the depth of the He II break at 4 
Ryd we can calculate the effective He II column density which 
gives rise to this break (Table 4). It is seen that the values 
of 7V(He II) are spread by maximum 40% around the median 
1.2 x 10 18 cm~ 2 . Besides, the recovered ionizing spectra show 
a diversity of shapes with a tendency to high spectral indices at 
1 < E < 4 Ryd (a > -1.4, median a = -1.0). These spectra 
are likely to be intrinsic spectra of quasars, i.e. not processed 
by the IGM. In order to cause a noticeable absorption in the 
He II Lyman continuum, the intervening absorbers should have 
iV(Hell) > 5 x 10 17 cm~ 2 . According to Fig. 19, the corre- 
sponding lower limit on N(H I) will be 3.0 x 10 15 cm -2 for the 
HM spectrum at z = 3 (middle curve) and 1.8 x 10 15 cm~ 2 and 



3 It is interesting to note that other absorbers from our study with 
characteristics similar to those of the complex at z = 2.735 — 
2.741, namely the 2 abs = 2.917 system towards HE 0940-1080 
(Sect. 3.1, Paper I) and the z abs = 2.875 system towards Q 1 157+3143 
(Sect. 3.2.3, present paper), also lie at the boundary of a void. 
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Fig. 18. Redshift dependence of the He II Lya opacity, tq^ 11 (2) . Open 
symbols show direct measurements from HE 2347-4342 (Zheng et al. 
2004), HS 1700+6416 (open square - Davidsen et al. 1996, crossed 
open square - Fechner et al. 2006a), and Q 1157+3143 (Reimers et al. 
2005b). Filled symbols represent the values of TQp 11 estimated from 
the restored spectral shapes of the UVB ionizing radiation: diamond 
and reverse triangle — Paper I; triangles, squares and circles — the 
present paper. Error bars correspond to 1 a confidence level. 
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Fig. 19. Dependence of 77 mo dci on the ionization parameter U for dif- 
ferent spectral shapes of the ionizing radiation: HM at z = 3 (solid 
line), and UVB restored at z — 2.93 (dotted line) and z — 1.80 
(dashed line) which are, respectively, softer and harder than the HM 
spectrum. The UVB at z = 2.93 (Q 1157+3143) is shown in Fig. 9, 
and at z = 1.80 (HS 0747+4259) — in Fig. 20. The adopted metallic- 
ityisZ = O.O1Z . 



5.5 x 10 15 cm -2 for the softer and harder spectra (upper and 
lower curves, respectively). From statistics of the Lya forest 
(e.g., Kim et al. 2002) it is well known that strong H I absorbers 
with AT(Hl) > 10 15 cm -2 become more rare with decreasing 
redshift. This results in less absorption at the He II ionization 
edge leading to a harder UVB. To produce iV(He II) > 5 x 10 17 
cm~ 2 , harder spectra, in turn, require more stronger H I ab- 
sorbers which are even more rare. Thus, by z ~ 1.8 a large part 
of the QSO/AGN emitted radiation remains unprocessed. 

In support to this consideration the following observation- 
based arguments can be given. Using a combination of several 
hard X-ray surveys, Ueda et al. (2003) investigated cosmologi- 
cal evolution of the AGN luminosity function. They found that 
the comoving spatial density of most luminous AGNs (quasars) 
peaks at z ~ 2 and after that decreases rapidly, but spatial 
density of less luminous but orders of magnitude more numer- 
ous AGNs continues to rise till z < 1. On the other hand, a 
FUSE survey of the AGNs (z < 1) by Scott et al. (2004) shows 
that lower luminosity AGNs tend to have harder spectral slopes 
(median a = -0.6) in the range 650-1000 A (0.9-1.4 Ryd). 

Thus, we can conclude that a significant hardening of the 
UVB spectrum occurs between z = 2.38 and z — 1.80. 
According to Fig. 19, this hardening would result in decrease 
of -q. The redshift evolution of 77 reported in Fechner & Reimers 
(2006) allows to suggest that hardening starts even earlier, 
somewhere at z ~ 2.5. At the moment we cannot trace the 
UVB evolution between z = 2.5 and z = 1.80 because of 
scarcity of our absorbers in this redshift interval. This will be a 
topic of our future studies. 



5. Conclusions 

1 . Using the MCISS method based on realistic models of line 
formation in stochastic media and adequate computational 
technique, we restored the spectral shapes of the underly- 
ing ionizing radiation from optically thin metal-line sys- 
tems. The UVB spectral shape does fluctuate in the redshift 
range studied, 1.8 < z < 2.94. However, the fluctuations at 
z < 1.8 and 2.4 < z < 2.94 are ruled by different physical 
reasons. At 2.4 < z < 2.94, the fluctuations arise due to 
radiative transfer processes in the clumpy IGM, whereas at 
z < 1.8 the IGM becomes almost transparent both in the 
H I and He II Lyman continua and the variability of spec- 
tral shapes comes from diversity of spectral indices of the 
intrinsic QSO/AGN spectra. 

2. At 2.4 < z < 2.94, the recovered spectral shapes show 
the intensity depression between 3 and 4 Ryd which is 
interpreted as a manifestation of the He II Gunn-Peterson 
effect. The values of the mean Hell Lya opacity calcu- 
lated from the depth of the recovered GP trough corre- 
spond well (within 1-2<t) to Tq£ h directly measured from 
the H i/He II Lya forest towards the quasars studied. This 
concordance confirms the applicability of the MCISS pro- 
cedure for studying He II absorption in the IGM. 

3. The scatter of rj = TV (He Il)/7V(H I) and anti-correlation be- 
tween A^(H I) and 77 can be explained by the combined ac- 
tion of the spectral shape variability due to radiative transfer 
in the clumpy IGM and differences in the mean gas densi- 
ties in the absorbing clouds. Progressive hardening of the 
recovered UVB towards lower z explains the correspond- 
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1. 1. Agafonova et al.: UVB and He II absorption at 1.8 < z < 3.0 



Table 4. He II column densities estimated from the ionizing spectra 
shown in Fig. 20 




10 E, Ryd 



Fig. 20. The difference between the UVB predicted by HM at 2 = 1.8 
(thick line) and those restored by the MCISS procedure (thin lines). 
The data are taken from (a) — the present paper, quasars HE 2347- 
4342 and HS 1700+6416, (b) — Reimers et al. (2005a), quasar HE 
0141-3932, (c) — Paper II, quasar HS 0747+4259, and (d) — Haardt 
&Madau(1996). 

ing decrease of the median //-values reported in F06 and 
Z04. 

4. All recovered UVBs are produced by QSO/AGN radiation 
processed by the clumpy IGM and do not show features 
which can be attributed to the input of radiation from soft 
sources like starburst galaxies. The metal systems consid- 
ered in the present paper and in Paper I have overdensi- 
ties 10-100 implying that radiative transfer effects both in 
the cloud and in the IGM can be significant. This means 
that for regions with smaller overdensities the UVB should 
be harder in general. However, such regions (voids) can 
be bounded by filament structures optically thin in H I but 
thick in He II. In this case, the incident radiation is effec- 
tively shielded and the UVB in voids becomes softer. 
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